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We consider an extension of the finite-temperature Dyson-Schwinger equation 
(DSE) approach to heavy mesons and quarkonia and apply it to calculate the 

cross section for the 3/tjj breakup reaction J/rp + tt > D + D . We study the 

effects of chiral symmetry restoration in the light quark sector on this process and 
obtain a critical enhancement of the reaction rate at the chiral/ deconfinement 
transition. Implications for the kinetics of charmonium production in ultrarela- 
tivistic heavy-ion collisions are discussed with particular emphasis on the recently 
observed anomalous J/i/> suppression as a possible signal for quark matter forma- 
tion. 



1 Introduction 

Recent results of the NA50 collaboration at GERN SPS on anomalous J/i/j 
suppression in Pb-Pb collisions at 158 A GeVEJa have renewed the quest for 
a proper description of charmonium production in heavy-ion collisions. It has 
been emphasized that a threshold effect in the dependence of J/V> supres- 
sion like throne observed by NA50 could be a signal for quark-gluon plasma 
formation B'u'EI, see Fig. |l|. However, up to now there is no theory for the 
hadroproduction of charmonium and the present approaches have to be con- 
sidered as parametrizations of the data. We provide here a reparametrization 
of this threshold effect where we extract effective comover cross sections from 
the difference of the observed data points to the extrapolated nuclear absorp- 
tion within a generalized Glauber model 0, see Fig. ^, lower panel. This 
comover cross section shows an onset at the transverse energy of Et ~ 40 
GeV which we interpret as a signal for a critical phenomenon in dense mat- 
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Figure 1: Upper panel: 3/tp production in Pb-Pb collisions at 158 AGevii compared to 
the extrapolation from p-A and S-U collisions using a Glauber model (dashed line) with 
nuclear J /if) absorption. Lower panel: Comover absorption cross sections extracted from the 
difference of observed data to the nuclear absorption extrapolation. 



ter due to the chiral/deconfincmcnt transition. Within the present paper, we 
want to investigate the charmonium dissociation cross section within a DSE 
model of QCD at high temperature and densityQ and thus go beyond previous 
estimates of this quantity within a nonrelativistic potential model El. This is 
particularly important since pions as quasi Goldstone bosons are involved in 
the dissociation process and a nonrelativistic description of these states fails 
to describe their wave functions accurately. The DSE approach will be applied 
for the description of the temperature dependence of heavy meson masses, 
which can lead to a dramatic increase in the ^ dissociation cross section, as 
has been demonstrated previously Q. This behaviour could contribute to an 
interpretation of the still puzzling experimental situation. 

2 Dyson-Schwinger equation approach to light and heavy mesons 

The nonperturbative features of low energy QCD such as dynamical chiral sym- 
metry breaking and confinement can be successfully modeled by the truncation 
of the full DSE for the, renormalized dressed quark propagator S(p), see Fig. 0, 
to the rainbow levelQfl. This approach has been systematically extended to the 
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Figure 2: Self-energy diagram defining the fully dressed quark two-point function S(p). 

case of finite temperatures and chemical potentials^ within the imaginary 
time (Matsubara) formalism, where it gives insight into the interplay between 
chiral symmetry restoration and deconfinement found in Lattice QCD simu- 
lations 113 and predicts the phase transition parameters in the quark matter 
phase diagram Ej. 

Masses and decay properties of light mesons (n, p) have been studied using 
the finiteptemperature Bethe-Salpeter equations (BSEs) in the corresponding 
channels E2fEJ. The DSE approach has also proven to be successful for the 
description of heavy meson observables uJ which play a crucial role in the 
present work. 

The approach we present in the following will employ a confining, separa- 
odellij recently extended to finite temperatures in the low-mass sector 
In a rank-1 truncation the model gluon propagator reads D(p, q) = 
Dq ip{p 2 )ip(q 2 ), and the renormalized quark propagator takes the form 

S(p) = -ij^p^avip 2 ) + &s(p 2 ) (1) 

with av(p 2 ) = [p 2 + m 2 (p 2 )]^ 1 , as(p 2 ) = m (p 2 )&v(p 2 ) and the dynamical 
quark mass function is given by 

m{p 2 ) = m + Am(T, p)(p(p 2 ). (2) 

The solution of the quark DSE reduces to that of the gap equation 

Am(7» = ^°| J0- Mq 2 )(Ts[q 2 ) , ( 3 ) 

where a nonvanishing mass gap Am(T, fi) in the chiral limit m = signals 
spontaneous chiral symmetry breaking in the T, \i -plane of the phase diagram. 
We choose a Gaussian formfactor f(p 2 ) = exp(— p 2 /A 2 ) and obtain a fit of 
7r- and p- meson properties with Dq — 227.8 GeV~ 2 , A = 0.60 GeV, mo = 
7.2 MeV. The temperature dependence of the solution of the mass gap equation 
(H) results in a chiral restoration transition temperature T c = 150 MeV, see 
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upper panel of Fig. ^. The model is confining for all Am(T, p) > A/y/2e since 
it has then no solution of p 2 + m 2 (p 2 ) = for real p 2 (no quasiparticle poles). 
For the heavy quark propagators the limit m{p 2 ) = ijiq = const will be used 
where toq is the current mass of the heavy quark Ej. 

The BSE in ladder approximation for the mesonic bound states has the 
form 

X(P 2 )T M (q,P) = ~ J ^0-D[ q ,p)^S qi (j> + )T M (jp,P)S q2 (p-) ltl (4) 

with p + = p + -qP, p- = p — (1 — r])P, where 77 = 1/2 for m qi = m q2 and 
j] = 1 for m qi 3> m q2 . Here we will use the one-covariant meson vertex 
functions Tm — 7m -Vm ifM (p 2 ) , where for the pseudoscalar mesons 7^ = jjj — 
Id = *75 anc i f° r the vector mesons 7 P = ^d* = 7jm = 7m with the proper 
normalization constants Nm- For the present study we will use the exploratory 
ansatz (pm(p 2 ) — fip 2 )- The condition \(P 2 = ~M\ t ) = 1 leads to the meson 
mass formulae 

1 = l Do I (0?^ 2(p2) [v^&pHvv + v^&pKvs 



,(5) 



where the functions V^' S ^ (p, P) can be found elsewhereEl. With toq = m c = 
1.844 GeV the resulting masses for the heavy mesons (experimental value in 
brackets) are M D = 1.869 GeV ((1.8693 ± 0.0005) GeV), M D . = 2.006 GeV 
((2.0067 ± 0.0005) GeV) and M J/4> = 3.459 GeV ((3.09688 ± 0.00004) GeV). 



3 Triangle diagram for the D* — > Dir decay 

We consider this vector-pseudoscalar-pseudoscalar decay as a test for our rel- 
ativistic quark model. In the light meson sector this process describes the 
Ptttt decay which is an important ingredient of recent studies of dilepton 
production in heavy-ion collisions due to vector meson dominance—It has also 
been studied in the present model at finite temperature recently Il3. Here we 
calculate the D*Dir decay process as a quark loop diagram, see Fig. § The 
decay width of this process is given by 

Td * ^ ^ iq^T/s \\Ml„Ml,Ml) , (6) 



with the kinematic factor A(s, Mf, Mf ) = y/[s - (Mi + M 2 ) 2 ][s - (Mi - M 2 ) 2 \ 
The D*Dir coupling constant go* Dtt — \{A—B) is given by the functions A, B, 
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Figure 3: The quark loop diagram for D* 



Dn decay vertex. 



defined from the transition amplitude of this process 
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N c J J^Hp^d {[!,] [-t(pi • iWiA) + <rs(M [nd 
[~i(p2 -l)vv(pl) + °s(pl)] [ns] [~i(P3 • tWCp!) + ^s(p§)] }( 7 ) 



where T{p 2 ) = N D .N V N D (p 3 (p 2 ) and pi = p - P*/2, p 2 = V + P*/2, p 3 = 
P + -P71-/2 — Pd- Our result for the coupling constant gu-D-K = 10.7 agrees well 
with the experimental value g^F D „ = 10.0 ± 1.3 EE The next type of diagram 
to be calculated is the box diagram that occurs at fourth order in the meson 
field expansion. 

4 J/ip dissociation cross section as a quark exchange process 

The J/?/ 1 dissociation cross section in a hot and dense medium of correlated 
quarks (mesons) has been studied within a nonrelativistic potential modella for 
quark exchange (string- flip) processes. This calculation gives large energy de- 
pendent cross sections with a threshold enhancement (peak value about 6 mb) 
and an exponential tail. It is, however, questionable whether the inadequate 
treatment of the pions in these models spoils the results. The pions which are 
Goldstone bosons of the broken chiral symmetry should be described within 
a chiral Lagrangian model. The corresponding reformulation of the charmo- 
nium dissociation process given previously E2I has recently been extended and 
improved EiJ. These approaches disregard the quark substructure effects which 
give rise to the nonlocality of the formfactors and which are expected to be- 
come apparent at finite temperatures and densities close to the QCD phase 
transition. 

Therefore, we consider here the formulation of the J/ip dissociation process 
within a relativistic quark model at finite temperature. As a generic new 
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Figure 4: The (anomalous) J/ip dissociation process in a relativistic quark model. Quark 
exchange (box diagram) is found to dominate over D* exchange (ladder resummed box 
diagram). 



mechanism on the quark level that we will study in detail is the anomalous 
process J/ip + n — ► D + D is described by the two types of diagram shown in 
Fig. ||. One type represents the quark exchange process (box diagram) and the 
other one describes the D* exchange (ladder resummed box diagram). In the 
language of the effective meson Lagrangian model, the box diagram is a contact 
term which is found to dominate over the heavy-meson exchange. Therefore we 
will focus on the evaluation of the cross section for the anomalous box diagram 
as a function of the center of mass energy y/s — y/JPj m + P-k) 2 



a(s) 



dt 



iTj/i/in dd\ 

16n\(s,M^M 2 ) 



(8) 



with the transition amplitude for this process being T 



J/ip -kD D 



and the kine- 



matic function A as defined above. In our calculations we put = 0, so 
that A = [s — Mj,^], s > 4M|, , and the physical region for this process is 
determined by 
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(9) 



In the DSE quark model the transition amplitude is given by 

d A p 1 



T 



J/t/jn DD 



Xtrr,|[-i(pi • 7)cry + crs}^D [~«(P2 • 7) + m Q] T j/4, 



-i{P3 ■ 7) + mq] T D [-i(p 4 ■ j)a v + (j s ] 



(10) 
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where Pl = p + P 7T /2,p 2 = p + P„/2 - P D ,p a = p- P n /2 + P D ,P4 = P - P-n/2. 
Our result for the s dependence of the cross section cr^™^ D5 is shown in Fig. 
§. We see a critical enhancement just above the reaction threshold. Then the 
cross section drops down for higher values of s. This qualitative behavior ia 
similar to the result of a calculation within a non-relativistic potential modelu 
but the absolute value of the cross section without in-medium effects (n = 0) 
is much smaller since we have considered the anomalous process first which 
vanishes in the non-relativistic treatment. 




Figure 5: Dependence of the J/1/1 dissociation cross section for the anomalous process (box 
diagram) on the squared center-of-mass energy s and on the parton density n of the medium. 

Now we investigate the influence of a hot and dense medium on char- 
monium dissociation. A first effect is the modification of the D-meson mass 
Mjj(T) due to chiral symmetry restoration. Similar to the temperature de- 
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pendence of the light quark mass the D-mass is nearly constant for lower 
temperatures and becomes smaller in the vicinity of the critical temperature, 
see Fig. |(| Secondly, this effect of chiral symmetry restoration enhances the 
cross section for the breakup reaction J/tjj + tt — > D + D. The enhancement 
of the breakup cross section has a critical temperature/density dependence at 
the chiral restoration transition. 



5 J/tjj kinetics at chiral symmetry restoration 

Aiming at a qualitative discussion of possible observable effects of the chiral 
restoration transition to be seen in the 3/tp production cross section in ultrarel- 
ativistic heavy-ion collisions we consider a relaxation time approximation for 
the evolution of the 3/tp distribution in a hot and dense pion gas with the pion 
distribution function f^(p;T,fi v ) = {exp[(^/|pl 2 + M% - ^)/T] - l}" 1 . This 
distribution function increases strongly near and above the critical tempera- 
ture. This leads to a critical enhancement of thermally averaged charmonium 
dissociation cross section (crj/^ i> rc i) and to a critical drop of the relax- 
ation time r 



T 



of the J/ij) distribution in a comoving dense pion (parton) gas. The result is 
shown in Fig. ^ lower panel. The critical enhancement of the thermally aver- 
aged dissociation cross section at the QCD phase transition temperature may 
lead to the result that in a heavy-ion collision above a critical energy density 
(Et) threshold this additional absorption process is switched on and results in 
an enhanced J /ip suppression. This pattern may serve as an explanation for the 
puzzling observation of anomalous 3/ip suppression in the NA50 experiment 
at CERN SpS, see Fig. 0. 



6 Conclusions 

A recently developed dynamical approach to the chiral symmetry restoration 
and deconfincmcnt transition is generalized for heavy mesons and quarkonia; 
in-medium modifications of masses and decay constants are studied. We apply 
this model to calculate the cross section for charmonium dissociation by hadron 
impact and find that contributions from the quark exchange dominate over 
those from D* resonance propagation processes. At the QCD phase transition 
the reaction rates for the 3/ip breakup are strongly enhanced. 

The approach will be further developed in order to include the p-meson 
which catalyzes exothermic 3/tp breakup and a treatment of the normal pro- 
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Figure 6: Upper panel: Double quark mass gap (solid line) and D meson mass (dashed line) 
as a function of temperature. Lower panel: Temperature dependence of the average cross 
section (arbitrary units) for J/ip dissociation (anomalous process) in a hot pion gas with 
critical enhancement at the chiral restoration temperature T c = 150 MeV. The dashed band 
shows the influence of a finite (nonequilibrium) pion chemical potential. 

cess J/if> + 7r — > DD* , D* D, which should be larger than the anomalous one 
considered here. On this basis, a detailed quantitative calculation of 3/tp pro- 
duction in the NA50 experiment can be attacked and predictions for the RHIC 
and LHC colliders can be made. A particularly new aspect is the occurrence of 
DD annihilation into J/ip + ir(p) CgaixL process) under these conditions which 
modifies the charmonium kinetics E3Ej and may lead to a saturation or even 
enhancement of the J ftp abundance. 
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